Sudden vanishing of spin squeezing under decoherence 
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In order to witness multipartite correlations beyond pairwise entanglement, spin-squeezing param- 
eters are analytically calculated for a spin ensemble in a collective initial state under three different 
decoherence channels. It is shown that, in analogy to pairwise entanglement, the spin squeezing de- 
scribed by different parameters can suddenly become zero at different vanishing times. This finding 
shows the general occurrence of sudden vanishing phenomena of quantum correlations in many- 
body systems, which here is referred to as spin-squeezing sudden death (SSSD). It is shown that 
the SSSD usually occurs due to decoherence and that SSSD never occurs for some initial states in 
the amplitude-damping channel. We also analytically obtain the vanishing times of spin squeezing. 
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I. INTRODUCTION 

Quantum entanglement [l| plays an important role in 
both the foundations of quantum physics and quantum- 
information processing [2]. Moreover, various entangled 
states have been produced in many experiments for dif- 
ferent goals when studying various nonclassical phenom- 
ena and their applications (3l-[l]|. Thus, entanglement 
is a quantum resource, and how to measure and detect 
entanglement is very crucial for both theoretical investi- 
gations and potential practical applications. 

For a system of two spin- 1/2 particles or a compos- 
ite system of a spin- 1/2 and a spin-1, there are opera- 
tionally c omp utable entanglement measures such as con- 
currence dl] and negativity [H, [l3|, but no universal 
measures have been found for general many-body sys- 
tems. To overcome this difficulty, entanglement witnesses 
are presented to detect some kinds of entanglement in 
many-bod y sy stems (T^ . [T5| . Now it is believed that spin 
squeezing [la, [l^ niay be useful for this task [l8l - [20j . In 
a general sense, spin-squeezing parameters are multipar- 
tite entanglement witnesses. For a class of many-particle 
states, it has been proved that the concurrence is linearly 
related to some squee zing pa rameters ^21i |. In fact, spin- 
squeezing parameters |16l - [l9| could be calculated also in 
a simple operational fashion, which characterizes multi- 
partite quantum correlations beyond the pairwise entan- 
glement. Another important reason for choosing spin- 
squeezing parameters as indicators of multipartite corre- 
lations is that spin squeezing is relatively easy to gener- 
ate [U [12] and measure experimentally [1^, [23| • 

Besides being a parameter characterizing multipartite 
correlations, spin squeezing is physically natural for con- 
trolling many-body systems. It is difficult to control a 
quantum many-body system since its constituents can- 



not be individually addressed. In this sense, one needs 
to use collective operations, and spin squeezing is one of 
the most successful approaches for controlling such sys- 
tems. For example, creating spin squeezing of an atomic 
ensemble could result in precision measurements based 
on many-atom spectroscopy 17[. Therefore, we can also 



regard spin squeezing as a quantum resource since for 
more than two particles it behaves as two-particle en- 
tanglement in controlling and detecting quantum corre- 
lations. On this quantum resource, we need to further 
consider the effects of decoherence [1^, [1^. Thus, it is 
important to study the environment-induced decoherence 
effects on both spin squeezing and multipartite entan- 
glement [27I - I37I . A decaying time evolution of the spin 
squeezing under decoherence [13, 138l - l40| can be used to 
analyze whether this quantum resource is robust. 

In this article we address this problem by calculating 
three spin-squeezing parameters for a spin ensemble in a 
collective excited state. We study the time evolution of 
spin squeezing under local decoherence, acting indepen- 
dently and equally on each spin. Here, the irreversible 
processes are modelled as three decoherence channels: 
the amplitude damping, pure dephasing and depolariz- 
ing channels. We find thatj_similar to the sudden death 
of pairwise entanglement |4l| . spin squeezing can also 
suddenly vanish with different lifetimes for some deco- 
herence channels, showing in general different vanishing 
times in multipartite correlations in quantum many-body 
systems. Thus, similar to the discovery of pairwise en- 
tanglement sudden death (BSD) j41|], the spin-squeezing 
sudden death (SSSD) occurs due to decoherence. We will 
see that for some initial states, the SSSD never occurs un- 
der the amplitude-damping channel. We also give analyt- 
ical expressions for the vanishing time of spin squeezing 
and pairwise entanglement. The ESD has been tested ex- 
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perimentally [39|, |42| and we also expect that the SSSD 
can also be realized experimentally. 

This article is organized as follows. In Sec. II, we intro- 
duce the initial state from the one-axis twisting Hamil- 
tonian and then, in Sec. Ill, the decoherence channels. 
In Sec. IV, we list three parameters of spin squeezing 
and discuss the relations among them. For a necessary 
comparison, the concurrence is also calculated. We also 
study initial-state squeezing. In Sec. V, we study three 
different types of spin squeezing and concurrence under 
three different decoherence channels. Both analytical and 
numerical results are given. We conclude in Sec. VI. 



II. INITIAL STATE 

We consider an ensemble of N spin- 1/2 particles with 
ground state |1) and excited state |0). This system has 
exchange symmetry, and its dynamical properties can be 
described by the collective operators 



N 

J a — ^ ^ jkc 

k=l 



1 ^ 

2 ^ 



(1) 



for a = x,y,z. Here, aka are the Pauli matrices for the 
kth qubit. To study the decoherence of spin squeezing, 
we choose a state which is initially squeezed. One typical 
class of such spin-squeezed states is the one-axis twisting 
collective spin state [l6l |. 

\^^J{9o))o = e-''°'^"/^\l)^^ = e-'^°'^"/^\l), (2) 

which could be prepared by the one-axis twisting Hamil- 
tonian 



where 



00 = 2xt 



(3) 



(4) 



is the one-axis twist angle and y is the coupling constant. 
For this state, it was proved [2l| that the spin squeezing 
[ll] and the concurrence Co [ll] are equivalent since 
there exists a linear relation 



state. Since Ja define an angular-momentum spinor 
representation of SO (3), the general definitions of spin 
squeezing for abstract operators Jx, Jy, and Jz can work 
well by identifying N/2 with the highest weight J, which 
corresponds to the collective ground state 



|J,-J) = |1 



(5) 



indicating that all spins are in the ground state. The 
symmetric space is generated by the collective operator 

1 ^ 

fc=i 

acting on the collective ground state. Here, 

. In others words, the state is in the maximally symmet- 
ric space spanned by the Dicke states. So, the N spin-1/2 
system behaves like a larger spin-7V/2 system. It can be 
proved that any pure state with exchange symmetry be- 
longs to the above-mentioned symmetric space, but for 
mixed states the state space can be extended to include 
a space beyond the symmetric one In the following 
discussions, we focus on such an extended space. 

In fact, after decoherence, not only the symmetric 
Dicke states are populated, but also states with lower 
symmetry. So, it is not sufficient to describe the sys- 
tem in only [N + l)-dimensional space. Although the 
maximal symmetry is broken, the exchange symmetry is 
not affected by the decoherence as each local decoherence 
equally acts on each spin. In other words, a state with 
exchange symmetry does not necessarily belong to the 
maximally symmetric space. 

With only the exchange symmetry, from Eq. ([T|), the 
global expectations or correlations of collective operators 
are obtained as 



{Jl) 



N N{N-1) 
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(6) 
(7) 
(8) 



between them. Physically, they occur and disappear si- 
multaneously. The spin squeezing of this state can be 
generated and stored in, e.g., a two-component Bose- 
Einstein condensate 14311. 



A. Initial-state symmetry 



Furthermore, it follows from Eq. ([B]) that 



2 N , N{N~l) , 

•Jy) = 2 {<^l+(^2- 



N \ N 



CTi • CT2) 



(9) 

(10) 



The initial state has an obvious symmetry resulting 
from Eq. the so-called even-parity symmetry, which 
means that only even excitations of spins occur in the 



These equations show the relations between the global 
and local expectations and correlations, which are useful 
in the following calculations. 
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III. DECOHERENCE CHANNELS AND 
EXAMPLES OF THEIR IMPLEMENTATIONS 

Having introduced the initial state, now we dis- 
cuss three typical decoherence channels: the amplitude- 
damping channel (ADC), the phase-damping channel 
(PDC), and the depolarizing channel (DPC). 

These channels are prototype models of dissipation rel- 
evant in various experimental systems. They provide 
"a revealing caricature of decoherence in realistic phys- 
ical situations, with all inessential mathematical details 
stripped away" But yet this "caricature of deco- 

herence" leads to theoretical predictions being often in 
good agreement with experimental data. Examples in- 
clude multiphoton systems, ion traps, atomic ensembles, 
or a solid-state spin systems such as quantum dots or 
NV diamonds, where qubits are encoded in electron or 
nuclear spins. 

Here, we briefly describe only a few of such implemen- 
tations. 



A. Amplitude-damping channel 

The ADC is defined as 

fADc(p) = EopEl + EipEl (11) 

where 

i?o = V^|0)(0| + |1)(1|, i?i = Vp|l)(0| (12) 

are the Kraus operators, p = 1 — s, s — exp(— 7i/2), 
and 7 is the damping rate. In the Bloch representation, 
the ADC squeezes the Bloch sphere into an ellipsoid and 
shifts it toward the north pole. The radius in the xy 
plane is reduced by a factor ^/s, while in the z direction 
it is reduced by a factor s. 

The ADC is a prototype model of a dissipative interac- 
tion between a qubit and its environment. For example, 
the ADC model can be applied to describe the sponta- 
neous emission of a photon by a two-level system into an 
environment of photon or phonon modes at zero (or very 
low) temperature in (usually) the weak Born-Markov ap- 
proximation. The ADC can also describe processes con- 
tributing to Ti relaxation in spin resonance at zero tem- 
perature. Note that by introducing an "upward" decay 
(i.e, a decay toward the south pole of the Bloch sphere), 
in addition to the standard "downward" decay, the ADC 
can be used to describe dissipation into the environment 
also at finite temperature. 

The ADC acting on a system qubit in an unknown 
state p can be implemented in a two-qubit circuit per- 
forming a rotation Ry{9) of an ancilla qubit (initially in 
the ground state) controlled by the system qubit and fol- 
lowed by a controUed-NOT (CNOT) gate on the system 
qubit controlled by the ancilla qubit [2]. The parameter 
is simply related to the probability p in Eq. (|lip . The 



ancilla qubit, which models the environment, is measured 
after the gate operation. 

The ADC-induced sudden vanishing of entanglement 
was first experimentally demonstrated for polarization- 
encoded qubits |42i]. For this reason let us shortly de- 
scribe this optical implementation of the ADC. It is based 
on a Sagnac-type ring interferometer composed of a po- 
larizing beam splitter and a half-wave plate at an angle 
corresponding to the parameter p in Eq. (ITU) . The beam 
splitter separates an incident beam (being in a superpo- 
sition of states with horizontal, \H), and vertical, \V), 
polarizations) into spatially distinct counter propagating 
light beams. The H component leaves the interferome- 
ter unchanged. But the V component is rotated in the 
wave plate, which corresponds to probabilistic damping 
into the H component. Then, at the exit from the inter- 
ferometer, this component is probabilistically transmit- 
ted or reflected from the beam splitter. So it is cast into 
two orthogonal spatial modes corresponding the reservoir 
states with and without excitation. 

The action of the ADC can be represented by an inter- 
action Hamiltonian 0]: iJ ab'^ + a^b, where a (a^) and 
b {h^) are annihilation (creation) operators of the sys- 
tem and environment oscillators, respectively. In more 
general models of damping, a single oscillator b of the 
reservoir is replaced by a finite or infinite collection of 
oscillators coupled to the system oscillator with dif- 
ferent strengths (see, e.g., Ref. 0,|43|). For the example 
of quantum states of motion of ions trapped in a radio- 
frequency (Paul) trap, the amplitude damping can be 
modeled by coupling an ion to the motional amplitude 
reservoir described by the above multioscillator Hamil- 
tonian [4^. The high-temperature reservoir is possible 
to simulate by applying (on trap electrodes) a random 
uniform electric field with spectral amplitude at the ion 
motional frequency 48, 49] . The zero-temperature reser- 
voir can be simulated by laser cooling combined with 
spontaneous Raman scattering (50| . 



B. Phase-damping channel 

The PDC is a prototype model of dephasing or pure 
decoherence, i.e., loss of coherence of a two- level state 
without any loss of system's energy. The PDC is de- 
scribed by the map 

W(p) -sp + p(poo|0)(0|+pii|l)(l|), (13) 

and obviously the three Kraus operators are given by 

E^ = ^st, = v^|o)(o|, £;2 = v^|i)(i|, (14) 

where 1 is the identity operator. For the PDC, there is 
no energy change and a loss of decoherence occurs with 
probability p. As a result of the action of the PDC, the 
Bloch sphere is compressed by a factor (1 — 2p) in the xy 
plane. 
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In analogy to the ADC, the PDC can be considered as 
an interaction between two oscillators (modes) represent- 
ing system and environment as described by the interac- 
tion Hamiltonian: H ^ a^a{b^ +b) 0. In more general 
phase-damping models, a single environmental mode b 
is usually replaced by an infinite collection of modes 6„ 
coupled, with various strengths, to mode a. 

It is evident that the action of the PDC is nondissipa- 
tive. It means that, in the standard computational basis 
|0) and the diagonal elements of the density matrix 
p remain unchanged, while the off-diagonal elements are 
suppressed. Moreover, the qubit states |0) and |1) are 
also unchanged under the action of the PDC, although 
any superposition of them (i.e., any point in the Bloch 
sphere, except the poles) becomes entangled with the en- 
vironment. 

The PDC can be interpreted as elastic scattering be- 
tween a (two-level) system and a reservoir. It is also a 
model of coupling a system with a noisy environment via 
a quantum nondemolition (QND) interaction. Note that 
spin squeezing of atomic ensembles can be generated via 
QND measurements [H H E^^- So modeling the 
spin-squeezing decoherence via the PDC can be relevant 
in this context. 

The PDC is also a suitable model to describe T2 relax- 
ation in spin resonance. This in contrast to modeling Ti 
relaxation via the ADC. 

A circuit modeling the PDC can be realized as a simpli- 
fied version of the circuit for the ADC, discussed in the 
previous subsection, obtained by removing the CNOT 
gate pj. Then, the angle 6 in the controlled rotation 
gate Ry{0) is related to the probability p in Eq. (fT3)) . 

The sudden vanishing of entanglement under the PDC 
was first experimentally observed in Ref. (43 |. This op- 
tical implementation of the PDC was based on the same 
system as the above-mentioned Sagnac interferometer for 
the ADC but with an additional half- wave plate at a 7r/4 
angle in one of the outgoing modes. 

Some specific kinds of PDCs can be realized in a more 
straightforward manner. For example, in experiments 
with trapped ions, the motional PDC can be imple- 
mented just by modulating the trap frequency, which 
changes the phase of the harmonic motion of ions [i^, 
(for a review see Ref. [i^l and references therein) . 



C. Depolarizing channel 



The definition of the DPC is given via the map 



where 



fDPc(p)=^i?fcp£^l, 



(15) 



j=0 




(16) 



are the Kraus operators. By using the following identity 



c^xP Ox + Oyp Gy + a^paz 



21, 



we obtain 



'^^DPc(/o) ^ sp + p 



1 



(17) 



= (1 - p')p + —{(J^pa^ + (Jypay + azpdz), , 



where p = 4p'/3. We see that for the DPC, the spin is 
unchanged with probability s = 1 — p or it is depolarized 
to the maximally mixed state 1/2 with probability p. It 
is seen that due to the action of the DPC, the radius of 
the Bloch sphere is reduced by a factor s, but its shape 
remains unchanged. 

Formally, the action of the DPC on a qubit in an un- 
known state p can be implemented in a three-qubit circuit 
composed of two CNOT gates with two auxiliary qubits 
initially in mixed states 

Pi = l/2, p2 = (l-p)|00)(00|+p|ll)(ll|, (18) 

which model the environment. Qubit p2 controls the 
other qubits via the CNOT gates [2| . 

The DPC map can also be implemented by applying 
each of the Pauli operators [1, ax, (Jy, Uz] at random with 
the same probability. Using this approach, optical DPCs 
have been realized experimentally both in free space [s^ 
and in fibers [57| . where qubits are associated with po- 
larization states of single photons. In Ref. [s^, the DPC 
was implemented by using a pair of equal electro-optical 
Pockels cells. One of them was performing a ax gate and 
the other a dy gate. The simultaneous action of both ax 
and ay corresponds to a ay gate. The cells were driven 
(with a mutual delay of t/2) by a continuous-wave peri- 
odic square-wave electric field with a variable pulse du- 
ration r, so the total depolarizing process lasted 2r for 
each period. 

Analogous procedures can be implemented in other 
systems, including collective spin states of atomic ensem- 
bles. The coherent manipulation of atomic spin states by 
applying off-resonantly coherent pulses of light is a ba- 
sic operation used in many applications [58j]. We must 
admit that the standard methods enable rotations in the 
Bloch sphere of only classical spin states (i.e., coherent 
spin states). Nevertheless, recently ^2^ an experimental 
method has been developed to rotate also spin-squeezed 
states. 

It is worth noting that in experimental realizations of 
decoherence channels (e.g, in ion-trap systems [111), suf- 
ficient resources for complete quantum tomography are 
provided even for imperfect preparation of input states 
and the imperfect measurements of output states from 
the channels. 
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IV. SPIN-SQUEEZING DEFINITIONS AND 
CONCURRENCE 



Now, we discuss several parameters of spin squeezing 
and give several relations among them. To compare spin 
squeezing with pairwise entanglement, we also give the 
definition of concurrence. We notice that most previous 
investigations on ESD of concurrence were only carried 
out for two-particle system rather than for two-particle 
subsystem embedded in a larger system. For the ini- 
tial states, spin-squeezing parameters and concurrence 
are also given below. 



A. Spin-squeezing parameters and their relations 

1. Definitions of spin squeezing 

There are several spin-squeezing parameters, but we 
list only three typical and related ones as follows 



2 _ 4(AJsJ 

min 



N 



4(J)2 



(19) 
(20) 
(21) 



Here, the minimization in the first equation is over all 
directions denoted by n_L , perpendicular to the mean spin 
direction (J)/(J^)\ Xmin is the minimum eigenvalue of the 
matrix [l9| 



r = (7v-i)7 + c, 



(22) 



where 



lki^Cki-{Jk){Ji) for k,le{x,y,z} = {1,2,3}, (23) 



is the covariance matrix and C = [Cki] with 



Ckl = -^{JlJk + JkJl) 



(24) 



is the global correlation matrix. The parameters 
£,i,£,2, and ^1 were defined by Kitagawa and Ueda [l6j , 
Wineland et al. [T3l, and Toth et al. jJi], respectively. If 
^2 < 1 (?3 < 1)' spii^ squeezing occurs, and we can safely 
say that the multipartite state is entangled [3 [3 • Al- 
though we cannot say that the squeezed state via the 
parameter ^\ is entangled, it is indeed closely related to 
quantum entanglement [21| . 



2. Squeezing parameters for states with panty 

We know from Sec. II. A that the initial state has an 
even parity and that the mean spin direction is along 



the z direction. During the transmission through all the 
three decoherence channels discussed here, the mean spin 
direction does not change. For states with a well-defined 
parity (even or odd), the spin-squeezing parameter 
was found to be ^] 



- {{J! + J'y) - \{J')\) 



(25) 



Then, the parameter ^| given by Eq. (l20l) becomes 

^2 - 77TT2 ~ o/ 7 \2 ■ l-^"-' 



For the third squeezing parameter (see Appendix A for 
the derivation), we have 



m^^ij^) - 2N- 



where 



(27) 



(28) 



Note that the first parameter S^f becomes a key ingredient 
for the latter two squeezing parameters (^| and ^|). 



3. 



Spin-squeezing parameters in terms of local expectations 



For later applications, we now express the squeezing 
parameters in terms of local expectations and correla- 
tions, and also examine the meaning of <;^, which will be 
clear by substituting Eqs. ([T]) and ^ into Eq. (^5]). 

^ = 1+C,, 

- l + (iV-l)((ai,a2,)-(ai,)(a2,)). (29) 

Thus, the parameter is simply related to the corre- 
lation Czz along the z direction. A negative correlation 
Czz < is equivalent to < 1. It is already known that 
the spin-squeezing parameter S^f can be written as [60j 



= 1 + (^ - l)CfJ,H, 



(30) 



where C^^s^ is the correlation function in the direction 
perpendicular to the mean spin direction. So, the spin 
squeezing < 1 is equivalent to the negative pairwise 
correlations Ch^s^ < [60l |. 

Thus, from the above analysis, spin squeezing and neg- 
ative correlations are closely connected to each other. 
The parameter < 1 indicates that spin squeezing oc- 
curs along the z direction, and < 1 implies spin squeez- 
ing along the direction perpendicular to the mean spin 
direction. Furthermore, from Eq. (|27|) . a competition 
between the transverse and longitudinal correlations is 
evident. 
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By substituting Eqs. (O and ^ to Eq. ([25]), one can 
obtain the expression of in terms of local correlations 
(cri+CT2-) and (cri-(T2-) as follows: 

= l + (iV-l)(CTi + CT2- +cri-0-2+) 

-2(iV- l)|(ai_a2-)| 
= l + 2(7V-l)(ai+a2-)-|(ai_a2-)|). (31) 

The second equality in Eq. pip results from the exchange 
symmetry. From Eqs. ([1]), ([TU|) . and (f^. one finds 



^3^ = 



rin{e?,l+C,,} 



(1-7V-i)(cti -aa) +^-1' 



(32) 
(33) 



Thus, we have reexpressed the squeezing parameters in 
terms of local correlations and expectations. 



4-. New spin-squeezing parameters 

In order to characterize spin squeezing more conve- 
niently, we define the following squeezing parameters: 



Cfc = max(0, 



1 



a) 



k £ {1,2,3}. 



(34) 



This definition is similar to the expression of the concur- 
rence given below. Spin squeezing appears when > 0, 
and there is no squeezing when C| vanishes. Thus, the 
definition of the first parameter has a clear meaning, 
namely, it is the strength of the negative correlations as 
seen from Eq. ([50]) . The larger is Ci, the larger is the 
strength of the negative correlation, and the larger of is 
the squeezing. More explicitly, for the initial state, we 
have ^1 = I - {N ~ l)Co [2l|, so Ci is just the rescaled 
concurrence Cf = ^(0) = (A^ - l)Co M- 

Here, we give a few comments on the spin-squeezing 
parameter which represents a competition between 
and (ciz)^: the state is squeezed according to the 
definition of ^| if < (criz)^ We further note that ^ 



1 - 2El, 



(35) 



where is the linear entropy of one spin and it can be 
used to quantify the entanglement of pure states |14l |. So, 
there is a competition between the strength of negative 
correlations and the linear entropy 2El in the parameter 
and Ci > 2El implies the appearance of squeezing. 



calculated from the reduced density matrix. It is defined 
as [H 



C = max(0, Ai — A2 — A3 — A4 



(36) 



where the quantities A^ are the square roots of the eigen- 
values in descending order of the matrix product 



£'12 = Pl2(o-ly ® 0-2a)Pl2(0'la ® '^■2y)- 



(37) 



In (1571) . /9*2 denotes the complex conjugate of pi2- 

The two-spin reduced density matrix for a parity state 
with the exchange symmetry can be written in a block- 
diagonal form 63] 



P12 



w+ u 
u w_ 



w y 

y w 



(38) 



in the basis {|00), |11), |01), |10)}, where 



v± 



u 

y 



-(l±2(ai,) + (ai,(T2.)), (39) 



^ (1 - {<Tlz(T2z)) 

{(Jl+(J2+), 

(0-1 + 0-2-)- 



(40) 

(41) 
(42) 



The concurrence is then given by [6J] 

C = max{0,2(|u| - w) ,2(y- 7^^:^)} . (43) 

From the above expressions of the spin-squeezing pa- 
rameters and concurrence, we notice that if we know 
the expectation {(Tiz), and the correlations (0-1+0-2-), 
(o-i_o-2-)j and {(Tiz(72z), all the squeezing parameters and 
concurrence can be determined. Below, we will give ex- 
plicit analytical expressions for them subject to three de- 
coherence channels. 

C. Initial-state squeezing and concurrence 

We will now investigate initial spin squeezing and pair- 
wise entanglement by using our results for the spin- 
squeezing parameters and concurrence obtained in the 
last subsections. We find that the third squeezing pa- 
rameter ^1 is equal to the first one £,f. The squeezing 
parameter ^1 is given by (see Appendix B): 



^?(o)-i-a(o) 

= l-(7V-l)Co, 

= l-2{N-l){\uo\~yo), 



(44) 



B. Concurrence for pairwise entanglement 

It has been found that_the concurrence is closely re- 
lated to spin squeezing ;!2lj. Here, we consider its behav- 
ior under various decoherence channels. The concurrence 
quantifying the entanglement of a pair of spin- 1/2 can be 



where 

Co -i{[(l - COS^-2 Q^^2 ^ ^-^2 (g^/2) 

- 1 + cos^-^ 9o} 
is the concurrence 1211. 



.,27V-4 



{0o/2)Y- 
(45) 



7 



The parameter ^|(0) is easily obtained, as we know 
both ^1 (0) and (0-12)0 dSS]). For this state, following from 
Eq. (fTO|) . ((Ti • (72)0 — 1, and thus the third parameter 
given by Eq. ([55)1 becomes 

^32(0)=minK?(0),^2(0)] 

= min[{l-a(0),l+C22(0)], (46) 

where the correlation function is 

= ^ (1 + cos^"' ^o) - cos2^-2 (^?o/2) > 0. (47) 

The proof of the above inequality is given in Appendix 
C. 

As the correlation function Czz (0) and the concurrence 
Cr(0) are always > 0, Eq. (|46)) reduces to 

?3(0)=C?(0) = l-a(0). (48) 

So, for the initial state, the spin-squeezing parameters 
f|(0) and ^i(O) are equal or equivalently, we can write 
Ci (0) — CKO) = C'r(O) according to the definition of pa- 
rameter Q, given by Eq. ([34]). Below we made a summary 
of results of this section in Table I. 



V. SPIN SQUEEZING UNDER DECOHERENCE 

Now we begin to study spin squeezing under three dif- 
ferent decoherence channels. From the previous analysis, 
all the spin-squeezing parameters and the concurrence 
are determined by some correlation functions and expec- 
tations. So, if we know the evolution of them under deco- 
herence, the evolution of any squeezing parameters and 
pairwise entanglement can be calculated. 




0.1 0.2 0.3 0.1 0.2 0.3 

P V 



FIG. 1: (Color online) Spin-squeezing parameters ("I (red 
curve with squares), 1^3 (top green curve with circles), and 
the concurrence Cr (blue solid curve) versus the decoherence 
strength p = 1 — exp(— 7f) for the amplitude-damping chan- 
nel, where 7 is the damping rate. Here, So is the initial twist 
angle given by Eq. In all figures, we consider an ensem- 
ble of A'^ = 12 spins. Note that for a small initial twist angle 
^0 (e.g., So = O.Itt), the two squeezing parameters and the 
concurrence all concur. For larger values of So, the parame- 
ters (^2, Cli ^nd C become quite different and all vanish for 
sufficiently large values of the decoherence strength. 



B. Amplitude-damping channel 

1. Squeezing parameters 



A. Heisenberg approach 

We now use the Heisenberg picture to calculate the 
correlation functions and the relevant expectations. A 
decoherence channel with Kraus operators is defined 
via the map 

Then, an expectation value of the operator A can be cal- 
culated as (A) =Tt[A£{p)] . Alternatively, we can define 
the following map, 

£Hp)=J2kIpK,. (50) 

It is easy to check that 

{A)=TT[A£{p)]=Tr[£HA)p]. (51) 

So, one can calculate the expectation value via the above 
equation (j51l) . This is very similar to the standard 
Heisenberg picture. 



Based on the above approach and the Kraus operators 
for the ADC given by Eq. (|12p . we now find the evo- 
lutions of the following expectations under decoherence 



(see Appendix D for details) 

{(Jiz) = s{(Tiz)q -p, (52a) 

(cti-O-2-) = s((Ti_(T2-)o, (52b) 

(0-1+0-2-) = s(o-i+o-2-)o, (52c) 

{criz<J2z) = 5^(0-1^0-22)0 - 23p{aiz)a + p^. (52d) 



To determine the squeezing parameters and the concur- 
rence, it is convenient to know the correlation function 
Czz and the expectation (cti • (72), which can be deter- 
mined from the above expectations as follows: 

(Ol • (T2) =1 - .SJ5.T0, (53) 
Czz [{<yiz<J2z)o - (o-lz)o(o-2z)o) 

^s'CzziQ), (54) 

where 

= 1 + 2(0-2)0 + (0-120-22)0- (55) 
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TABLE I: Spin-squeezing parameters ^ [ij, ?! [3 and concurrence C ,12] for arbitrary states (first two columns), 

states with parity (third column). The squeezing parameters are also expressed in terms of local expectations (fourth column) 
and in terms of the initial rescaled concurrence Cr(0) for initial states (last column). Also, Co is the initial concurrence, and 
other parameters are defined in the text. 



Squeezing parameters 


Definitions 


States with parity 


In terms of local expectations 


Initial state 


SJ- 


4(AJfl^)Ln 

N 


— ((jj + jI) - \ {f)\) 


1 + 2(7V- l)(y- |u|) 


1 - Cr{0) 




2 






1 - a(o) 










S3 


-^min 

TV 




min{C?,l+C,4 


1 - a(o) 




(l-Af-i)(ai-a2) + iV-i 


Concurrence C 


max(0, Ai — A2 — A3 — A4) 


2 max(0, \u\ — w,y — yjv+v- ) 


2 max(0, \u\ — w,y — ^/v+v- ) 


Co 



Substituting the relevant expectation values and the cor- 
relation function into Eqs. (|3T|) . (|32p . and (pS)) leads to 
the exphcit expression of the spin-squeezing parameters 



ei = i-sa(o), 



nin{e^,l + g^C,,(0)} 
1 + {N-^ -l)spxo 



(56) 
(57) 

(58) 



As the correlation function Czz{0) > 0, given by Eq. (|47l) . 
the third parameter can be simplified as 



1 + -l)spxo' 



(59) 



Initially, the state is spin squeezed, i.e., ^i(O) < 1 or 
Cr{0) > 0. From Eq. one can find that ^1 < 1, 

except in the asymptotic limit of p = 1 . As we will see 
below, for the PDC and DPC, 



1 



'a(o). 



Thus, we conclude that according to the initially spin- 
squeezed state is always squeezed for p 7^ 1, irrespective 
of both the decoherence strength and decoherence mod- 
els. In other words, there exists no SSSD if we quantify 
spin squeezing by the first parameter ^f. 



2. Concurrence 

In the expression (j43|) of the concurrence, there are 
three terms inside the max function. The expression can 
be simplified to (see Appendix E for details): 



Cr = 2(7V- l)max(0, |u| - w). 



(60) 



By using Eqs. pO)) and (I52cp . one finds 

2{\u\~w) (61) 



2s|uo| + - [s-2 + s{aiza2z)o - 2^(0-1^)0]) 



spxo 



(62) 



So, we obtain the evolution of the rescaled concurrence 



Cr^max[0,sCriO)~2-\N -l)spxo] , (63) 

which depends on the initial concurrence, expectation 
{criz)o, and correlation {aiza2z)o- 



3. Numerical results 

The numerical results for the squeezing parameters and 
concurrence are shown in Fig. 1 for different initial values 
of the twist angle 9o, defined in Eq. For the smaller 
value of ^0, e.g., = tt/IO, we see that there is no ESD 
and SSSD. All the spin squeezing and the pairwise en- 
tanglement are completely robust against decoherence. 
Intuitively, the larger is the squeezing, the larger is the 
vanishing time. However, here, in contrast to this, no 
matter how small are the squeezing parameters and con- 
currence, they vanish only in the asymptotic limit. This 
results from the complex correlations in the initial state 
and the special characteristics of the ADC. 

For larger values of 6*0, as the decoherence strength p 
increases, the spin squeezing decreases until it suddenly 
vanishes, so the phenomenon of SSSD occurs. There ex- 
ists a critical value Pc, after which there is no spin squeez- 
ing. The vanishing time of ^| is always larger than those 
of ^1 and the concurrence. We note that depending on 
the initial state, the concurrence can vanish before or af- 
ter This means that in our model, the parameter 
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^1 < 1 implies the existence of pairwise entanglement, 
while ^1 does not. 



4- Decoherence strength pc corresponding to the SSSD 



From Eqs. ([57| . (p)) . and (|63 
can be analytically obtained as 



the critical value pc 



(k) 



XkCriO) 



,(2) 



iN~l)xo 

ha(o) 



(fc = 1,3) 



1 



l + 2((7l. 



\2 ' 

/o 



(64) 
(65) 



where xi = 2 for the concurrence and x^ = N for the 
squeezing parameter C|. The critical value pc^^ is for the 
second squeezing parameter C|. Here, Pc is related to the 
vanishing time tv via pc = 1 — exp(— 7it,). 

In Fig. 2, we plot the critical values Pc of the deco- 
herence strength versus ^o- The initial-state squeezing 
parameter is also plotted for comparison. For a range 
of small values of ^o, the entanglement and squeezing are 
robust to decoherence. The concurrence and parameter 
el intersect. However, we do not see the intersections 
between Q and C| or between Q and the concurrence. 
We also see that for the same degree of squeezing, the 
vanishing times are quite different, which implies that 
except for the spin-squeezing correlations, other type of 
correlations exist. For large enough initial twist angles 
TT < ^0 ^ 27r, the behavior of the squeezing parameter 
is similar to those corresponding to pi^^ and pi'^^ 



C. Phase-damping channel 

1. Squeezing parameters and concurrence 

Now, we study the spin squeezing and pairwise entan- 
glement under the PDC. For this channel, the expecta- 
tion values (erf") are unchanged and the two correlations 
(cri_cr2-) and ((Ti+(T2-) evolve as (see Appendix D for 
details) 



(cri_cr2-> = s^(CTi_cr2- 

((Tl + (T2-) = S^((Tl+(T2- 



(66) 



From the above equations and the fact (tJi • f72)o = 1, one 
finds 

(cti • f?2) — s^{p\x(yix + crij^cr2y)o + (o'l2Cr2z)o 

= S^(l - (CTi2Cr22)o) + (o-lzCr22)o, (67) 

CUp)=C^M- (68) 

Therefore, from the above properties, we obtain the evo- 
lution of the squeezing parameters. 



i-s^a(o). 



(69) 
(70) 




ADC 




0.5 1 1.5 

Initial twist angle Qoj'K 



FIG. 2: (Color online) Critical values of the decoherence 
strength p^P (blue solid curve), pP (red curve with squares), 
pP (green curve with circles), and the squeezing parameter 
("i (black dashed curve) versus the initial twist angle given 
by Eq. (U) for the amplitude-damping channel, ADC. Here, pc 
is related to the vanishing time t^ via pc = 1 — exp(— 7ti,). At 
vanishing times, SSSD occurs. The critical values pP , 
and p^c^ correspond to the concurrence, squeezing parameter 
Cl, and Cl, respectively. 



and the third parameter becomes 

A^min[^2^1-|-C^ 



3 (7v_i)[,,2 + (i_52)(^^^^2.)o] + l 

^ 

(Ar-l)[s2 + (l-s2)(^^^^2.)o] + l' 



(71) 
(72) 



where we have used Eqs. (|67|) and ((68)) . and the property 
C,M > 0. 

From Eq. (j66l) and the simplified form of the concur- 
rence given by Eq. (1601) . the concurrence is found to be 



max 



ix|o,2(7V- 1) 

[s^l^ol- 4-1(1 -(ai,(72.)o))] } 
ao(s2-l)" 



0,s^a(0) + 



where 



ao ^{N ((7i^a2.)o). 



(73) 



(74) 



Thus, we obtained all time evolutions of the spin- 
squeezing parameters and the concurrence. To study the 
phenomenon of SSSD, we below examine the vanishing 
times. 
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FIG. 4: (Color online) Same as in Fig. 2 but for the depolar- 
FIG. 3: (Color online) Same as in Fig. 2 but for the phase- j^ing channel, DPC, instead of ADC. 
damping channel, PDC, instead of ADC. 



Z. Decoherence strength pc correspond 



to the SSSD 



The critical decoherence strengths pc can be obtained 
from Eqs. dZO]), dZlI), and ([731) as follows: 



P^^' = 1 - 



XkCr{Q) + ao 



1 - (^i.>o 



a(o) 



(75) 
(76) 



where /c = 1, 3 and xi = 2, = TV. 

In Fig. 3, we plot the decoherence strength pc versus 
the twist angle 9q of the initial state for the PDC. For this 
decoherence channel, the critical value p'^s first decrease 
until they reach zero. Also, it is symmetric with respect 
to 00 = """j which is in contrast to the ADC. There are also 
intersections between the concurrence and parameter 
and the critical value p'i'^ is always larger than p^c^ and 



D. Depolarizing channel 

1. Squeezing parameters and concurrence 

The decoherence of the squeezing parameter defined 
by S0rensen et al. [11] has been studied in Ref. [l^l for 
the DPC. It is intimately related to the second squeezing 
parameter For the DPC, the evolution of correla- 
tions (cri_(T2-) and (cri+cr2-) are the same as those of 



the DPC given by Eq. and the expectations (ctiz) 

and {(Tiz<J2z) change as (see Appendix D). 



(0'lzCr2z 



5(0-12)0, 



S {Criz<72z)0- 

From these equations, we further have 



(ctI • CT2) = S^((Tl • (T2)o = S^, 



(77) 
(78) 



(79) 



Cz 



((ai.a2.)o - (fTi,)o(cT2z)o) = s^CzziO). (80) 



The squeezing parameter is given by Eq. (p^ . and the 
other two squeezing parameters are obtained as 



{N- l)s2 + 1 



{N ~ i)s2 + r 



(81) 



(82) 



By making use of Eqs. ( [521) and ([75]) and starting from 
the simplified form of the concurrence (|B(71) , we obtain 



Cr — max 
= max 



{0, 2{N - 1) [s^\uo\ - i(l - s^{<Ji.a2z)o)] } 
[0,s2c,(0) + 2-i(iV-l)(s2-i)] . (83) 



We observe that the concurrence is dependent only on 
the initial value itself, not other ones. 
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TABLE II: Analytical results for the time evolutions of all relevant expectations, correlations, spin-squeezing parameters, and 
concurrence, as well as the critical values of the decoherence strength p. This is done for the three decoherence channels 
considered in this work. For the concurrence C, we give the expression for C'r, which is related to the rescaled C0IlCU.rr6IlC6 Cr 
via Cr = max(0, C^). 





Amplitude-damping channel 
(ADC) 


Phase-damping channel 
(PDC) 


Depolarizing channel 
(DPC) 




s(o-lz)o — P 


(012)0 


3(012)0 




S^(o-lzO-2z)o - 2sp{(Ji^)o +p'^ 


(o-l20-2z)o 


5^(012022)0 


(0-1 + 0-2-) 


s(o-i+02_)o 


S^(oi + 0-2-)o 


S^(oi + 0-2-)o 


{0-I-O2-) 


s(oi_0-2-)o 


S^(o-i_0-2-)o 


S^(o-l_02-)o 


{Si ■ 02) 


I — spxo 


s^(l — (01^022)0) + (oi 
















1 - sCr{0) 


1 - S^Cr{0) 


1 - s^aio) 


& 


1 - sCr{0) 

(s(oiz)o ~pr 


1 - S^Cr{0) 
(oi,)2 


1 - S^Cr{0) 
s2(oi2)g 


?3 


1 - sCr(O) 


1 - S^Cr{0) 




1 - S^Cr{0) 




1 -1- (iV-1 - l)spxo 


(1 - Ar-i)[s2 + (1 - s2){ai,a2,)o] + N-^ 


(1 






sC,-(0) -{N- l)spxo/2 


s'^CriO) + ao{s'^ - l)/2 




L)/2 


Pc 


2a (0) 




^ ao ^ 


1 

' 


1 


( N-l \ 


1 
2 


{N-l)xo 


V2a(0) +ao y 


\2Cr{Q) + N - I ) 


(2) 
Vc 


+ CriO) - 1 


1 - 




1 

2 


1 - 




1 
2 


l + 2(criz)o + (oz)g 


V c.(o) ; 


(a(o) + (oi2)g) 


(3) 
Vc 


iva(o) 




ao 


)* 




A- 1 




(Af - 1) xo 




,A-a(o) + A - 1^ 



Decoherence strength pc corresponding to the SSSD 



From Eqs. ((83|) . (|8ip . and (j82|) . the vanishing times are 
analyticaUy calculated as 



pi'' 



N~ 1 



XkCr{0)+N -1 

1 

a(0) + 



(84) 
(85) 



■where fc = 1, 3 and xi = 2,X3 — N. 

In Fig. 3, -we plot the critical values Pc versus the initial 
t-wist angle 6*0 for the DPC. For the DPC, the p'^s first in- 
crease until they reach their maxima and then decrease 
to zero. Also, it is symmetric -with respect to = ""j 
-which is the same as for the PDC. There are also inter- 



sections bet-ween the concurrence and the parameter 

Qualitatively, the behaviors of pi^' and pi^' are the same 
as that of the squeezing parameter (^f . This implies that 
the larger the squeezing, the larger is the critical value 

Pc- 



The common features of these three decoherence chan- 
nels are: (i) The critical value p^s is al-ways larger or equal 
than the other t-wo, namely, the spin-squeezing correla- 
tions according to ^| are more robust; (ii) there al-ways 
exist t-wo intersections bet-ween the concurrence and the 
parameter for 9o from to 27r, irrespective of the 
decoherence channels; (iii) -when there is no squeezing 
(central area of Figs. 2, 3, and 4), all vanishing times are 
zero. Table II conveniently lists all the analytical results 
obtained in this section. 
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VI. CONCLUSIONS AND REMARKS 

To summarize, for a spin ensemble in a typical spin- 
squeezing initial state under three different decoherence 
channels, we have studied spin squeezing with three dif- 
ferent parameters in comparison with the pairwise entan- 
glement quantified by the concurrence. When the sub- 
systems of the correlated system decay asymptotically in 
time, the spin-squeezing parameter also decays asymp- 
totically in time for all three types of decoherence. How- 
ever, for the other two squeezing parameters C| and C|, 
we find the appearance of spin-squeezing sudden death 
and entanglement sudden death. The global behaviors of 
the correlated state are markedly different from the local 
ones. The spin-squeezing parameter C| can vanish be- 
fore, simultaneously, or after the concurrence, while the 
squeezing parameter C| is always the last to vanish. This 
means that this parameter is more robust to decoherence, 
and it can detect more entanglement than 

Our analytical approach for the vanishing times can be 
applied to any initial quantum correlated states, not re- 
stricted to the present one-axis twisted state. Moreover, 
for more complicated channels, such as the amplitude- 
damping channel at finite temperatures ^Slij or the chan- 
nel discussed in Ref. the method developed in this 
article can be readily applied to study spin squeezing un- 
der these decoherence channels. 

Our investigations show the widespread occurrence of 
sudden death phenomena in many-body quantum cor- 
relations. Since there exists different vanishing times 
for different squeezing parameters, spin squeezing of- 
fers a possible way to detect the total spin correlation 
and their quantum fluctuations with distinguishable time 
scales. The discovery of different lifetimes for various 
spin-squeezing parameters means that, in some time re- 
gion, there still exists another quantum correlation when 
other quantum correlations suddenly vanish. However, 
to determine which kind of correlations will vanish, one 
possible approach is to further invoke irreducible mul- 
tiparty correlations (66j . where the multipartite correla- 
tions are classified in a series of irreducible k party ones. 
If we could obtain the time evolution behaviors of such 
irreducible multipartite correlations in various decoher- 
ence channels, we could classify lifetimes for the spin- 
squeezing sudden death of various multipartite correla- 
tions order by order. 
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Appendix A: Spin-squeezing parameter for states 
with parity symmetry 



Here, we calculate the spin-squeezing parameter ^| for 
collective states with either even or odd parity symmetry. 
For such states, we immediately have 



{Jx) = {Jy) = {JxJz) = {JyJz) = 



(Al) 



'V/ ^ \'Jx'Jzl — X-Jy- 

as the operators change the parity of the state. Then, 
the mean spin direction is along the z direction and the 
correlation matrix given by Eq. (j24p is simplified to 



C = 



where C^y = {[Jx 




(A2) 



From the correlation ma- 



trix C and the definition of covariance matrix 7 given by 
Eq. (1121), one finds 



/ iV(jJ) NC^y 



r = 



V 







xy 



N{Jl) 









(A3) 



This matrix has a block-diagonal form and the eigenval- 
ues of the 2x2 block are obtained as 



Ah 



N 



{{Jl + Jl)±\{Jl)\) 



(A4) 



In deriving the above equation, we have used the relation 



Ji^ji-j'y-l[j,,jy\ + . (A5) 

mm of r is obtained 



Therefore, the smallest eigenvalue A 
as 



Aniin 



(A_,iV(AJ,)2 + (j2)) 



(A6) 



where A_ differs from the squeezing parameter ^\ given 
by Eq. (^51) by only a multiplicative constant, as seen 
by comparing Eqs. and (|A6|) . From Eqs. (|A6I) and 
(HH), one finds that the squeezing parameter ^| is given 
by Eq. 
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Here, we will use the Heisenberg picture to derive 
the relevant expectations and spin-squeezing parameters 
for the initial state 67, 68]. To determine the spin- 
squeezing parameter given by Eq. pip , one needs to 
know the expectation (1712)0, and correlations (cri+cr2-)o 
and ((Ti_cr2-)o- We first consider the expectation {a\z)o- 
For simplicity, we omit the subscript in the following 
formulas. 
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1. Expectation (aiz) 
The evolution operator can be written as, 

U = exp{-ixtJ^) = exp f -i9 ^ jkxjix j (Bl) 

\ k>l ) 

up to a trivial phase, where 9 = 2xt given by Eq. Q. 
From this form, the evolution of jiz can be obtained as 



U^JizU = ji, cos[0j(2)j ^ ^.^^ sm[9ji% (B2) 



where 



TV 



l=k 

Therefore, the expectations are 



(B3) 



(B4) 



since (l|jiy|l) = 0. Here, |1') = |1)2 ® ■■■ <E) \1)n- So, one 
can find the following form for the expectation values 

(l|cos[0J,]|l> - ((l|e'«^-|l)+c.c.)/2 

= (nf^i(l|e^«^-|l) + c.c.)/2 

= cos^(0'), (B5) 

where 9' ^ 9/2 and |1) = |1)®^. 

By using Eqs. (jB4l) and (IBSp . one gets 



(a, = - cos 



N-l fQl\ 



(B6) 



2. Correlation ((Ti+cr2-) 



Since the operator aixCf2x commutes with the unitary 
operator [/, we easily obtain 



(0'l2;O'2x) — 0. 



(B7) 



We now compute the correlations {aiz'J2z)- From the uni- 
tary operator, 

= ji,cos(^4^))+ji,sin(0j(2)) 

X [j2, cos[0(ji, + _^ -^^ sin[0(ji, + 43))] 

= ji2 cos{9j2x) cos{9ji^'>) - jiz sin{9j2x) sin{9jl^^^) 
+jiy sm{9j2x) cos{9ji^'>) + jiy cos{9j2x) sin(6'j(3)) 
X \]2z cos{9jix) cos(6'j(3)) - sin(6lji^ ) sin(6'j ) 
+32y sm{9jix) cos{9ji^'>) + j2y cos(0ji,) sin(0j(3)) 



Although there are 16 terms after expanding the above 
equation, only 4 terms survive when calculating {sizS2z)- 
We then have 

{jlzj2z) - (1|jW2.COs2(0/2)cOs2(^43)) 

-jizj2xj2y sm{9) sm'^{9ji^^ ) 

+4=jlyjlxj2xj2y siv?{9l2) C0s'^{9j^^^) 

-3iy3ixi2z siiv{9)siT?{9jf^)\\) 
= 4-1(1'|cos^(0jF))|1') 



l+COS(20jW) |;^,^ 



8"^ [1 + COS 



N~2f 



(B8) 



where = [1)3 ® ... (E) |l)Af- The second equality in 
Eq. (jBSP is due to the property jxjy = —jyjx = ijz/'^, 
and the last equality from Eq. (jBSp . Finally, from the 
above equation, one finds 



(fTl.fT2.) -2-1(1 +C0S^-2 



(B9) 



Due to the relation {(Tixa2x + '^iy<^2y + <^iz(^2z) = 1 for 
the initial state, the correlation {aiya2y) is obtained from 
Eqs. dlTl) and (|B9| as 



{<^ly(^2y) 



2-1(1 



.,N-2 



(BIO) 



Substituting Eqs. (|B7| and (jBlOp into the following re- 
lations 

<^lx(y2x + <^ly<y2y = 2 (o'i + cr2- + cri_0-2+) 

leads to one element of the two-spin reduced density ma- 
trix, 



2/0 



(0-1 + 0-2- 



-ifl-cos"-2( 



(Bll) 



where the relation ((Ti+(T2-) — (0-1-0-2+) is used due to 
the exchange symmetry. 



3. Correlation l^a\-G2-) 

To calculate the correlation (o-i_o-2-), due to the fol- 
lowing relations 

0-12:0-22: — O-iyO-21/ = 2 (^1 + 0-2+ + 0-i_ 0-2- ) , (B12) 
i {crix(^2y + O-iyO-22:) = 2 (0-1 + 0-2+ - 0-1-0-2-) , (B13) 

we need to know the expectations {jixj2y)- The evolution 
of jixj2y is given by 



U^SlxS2yU = jl2:|j2iyC0S 

- j2z sin 



0{jlx+ji'^)]}, 
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and the expectation is obtained as 
(jW2y) = 2-i(l'|ji,sin[0(ji,+j(3))] 
= (4*)-i(l'|ji,e^«^-nf^3e^''^- 

= (4i)-icos~-2(0')(l|ji,e'»J- -ji,e-^«^-|l) 



2-icos^-2 
4-isin (61')' 



(0')(lb-i.sin(0ji,)|l) 



Here, = |l)i (g) |1)3 ® ... ® |l)Ar, where |1)2 is ab- 
sent. Moreover, {jiyj2x) — {jixj2y) due to the exchange 
symmetry, and thus, 

(iixj2y + JiW2^) = 2-isin(0')cos^"' [0'). 

For the initial state ([2]), we obtain the following expecta- 
tions O, fH 



0'lyO'2x) 



2sin( 



(B14) 



The c ombination of Eqs. dUI, (|BT0| . (|BT2)) . (|BT31) . and 
(|B14[) leads to the correlation 



"0 = (o'i-Cr2-) = — 8 ^ (l — COS 

-i2-isin(6l')cos^-2 



Ar-2 , 



(B15) 



Substituting Eqs. ((BTT|) and (|BT5)) to Eq. dSj) leads to 
the expression of the squeezing parameter given by 
Eq. (glD. 



Appendix C: Proof of C^^(O) > 

To prove this, we will not use this specific function of 
the initial twist angle 9 as given by Eq. (|T7l) , but use the 
positivity of the reduced density matrix (I55|) . We first 
notice an identity 



which results from Eqs. ([39| and (|40|). This is a key step. 
Also there exists another identity 

Wo = yo (CI) 

as (iJi • 0^2)0 = 1- From the positivity of the reduced 
density matrix ([55]) . one has 

vo+vo- > |woP > 2/0 = "^0' 

where the second inequality follows from Eq. (l40l) and 
the last equality results from Eq. (jCl|) . This completes 
the proof. 



Appendix D: Derivation of the evolution of the 
correlations and expectations under decoherence 



For an arbitrary matrix 




from the Kraus operators ((T2)) for the ADC, it is straight- 
forward to find 



£{A) 



sa yfsb \ 
\/sc d + pa I 

sa + pd y/sb 



The above equations imply that 

SH^f^) = V^cr^j for pL = x.y, 



SCTz - P- 



As we considered independent and identical decoherence 
channels acting separately on each spin, the evolution 
correlations and expectations in Eqs. (j52bp . (|52cp . and 
(j52dp are obtained directly from the above equations. 

From the Kraus operators (fH|) . the evolution of the 
matrix A under the PDC is obtained as 

SiA) ^ 8\A) ^ (^^^ -^'^y 
from which one finds 

So expectations {erf") are unchanged and Eq. ([5^ is ob- 
tained. 

From the Kraus operators (fTB|) of the DPC, the evolu- 
tion of the matrix A is given by 

£iA) = S^A) 

(as + ^(a + d) sb \ 

sc ds + f (a + d) J 

from which one finds 

£\aa) — sUa for a € {x, y, z}. 
Then, Eq. ^ is obtained. 

Appendix E: Simplified form of the concurrence 

For our three kinds of decoherence channels, the con- 
currence (|43p can be simplified and given by 

C = max{0,2(|u|-u;),2(t/-/[^l^)} 

= max{0,2(|M| - w)} . (El) 



If one can prove 



H-y > 0, 

W — y/U^UZ < 0, 



(E2) 
(E3) 
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then we obtain the simplified form shown in Eq. (jEll) . 
The last inequality can be replaced by 

~ v+v^ < (E4) 

as w and are real. 

We first consider the ADC channel. From Eqs. ()52bl) . 
(|52cp . and (|54l) . one obtains 

\u\-y = s{\uo\ - yo) > 0, (E5) 
1 

- v+v^ = = -^C,,(0) < 0. (E6) 

where the inequalities result from Eqs. (j^J) and (|Tf|. 
respectively. So, the inequality (|E4p follows. 

For the PDC, from Eq. (jMl) and fact that (erf") is 
unchanged under decoherence, the concurrence can also 



be simplified due to the following properties: 

\u\-y = s^d^^ol - yo) > 0, 
w^-v+v- - -ic,,(0) <0. 

For the DPC, from Eqs. and (HE]), one has 

\u\-y = s^{\uo\-yo)>0, (E7) 
w^-v+v- = -— C,,(0)<0. (E8) 

So, again, the concurrence can be simplified to the form 
shown in Eq. (|Eip . This completes the proof. 
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